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Abstract: We calculate pion mass, pion decay constant, PCAC quark mass and nucleon mass 
in two flavour lattice QCD with unimproved Wilson fermion and gauge actions. Simulations are 
performed using DD-HMC algorithm at two lattice spacings and two volumes for several values 
of the quark mass. The cutoff effects in pion mass and nucleon mass for the explored region of 
parameter space are found to be negligible. The chiral behaviours of pion mass, pion decay constant 
and quark condensate are found to be qualitatively consistent with NLO chiral perturbation theory. 
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1 Introduction 

Because of the explicit violation of chiral symmetry by a dimension five kinetic operator, Wil- 
son formulation has been known to be difficult to simulate at light quark masses. Lack of chiral 
symmetry means that the "physical" quark mass is no longer proportional to the bare quark mass 
(the quark mass renormalization is no longer only multiplicative) and Wilson-Dirac operator is not 
protected from arbitrarily small eigenvalues and may lead to zero or near zero modes for individ- 
ual configurations. This is the infamous problem of "exceptional configurations". This leads to 
convergence difficulties for fermion matrix inversion. This poses difficulties for lattice simulations 
with Wilson fermions in the chiral region. 

In the past, simulations with unimproved Wilson action has shown large scaling violations 
in hardonic observables. However, one should keep in mind that most of these were quenched 
simulations done at large pion masses, not small enough lattice spacings and smaller volumes. 
Further, the demonstration of the suppression of topological susceptibility with decreasing quark 
mass was inconclusive. The chiral behaviour of pion mass and decay constant with respect to quark 
mass (specifically, the presence of chiral logarithms) as dictated by chiral perturbation theory was 
also not convincingly demonstrable in the past with un-improved Wilson fermions. All these issues 
raise the question: Does Wilson lattice QCD belong to the same universality class as continuum 
QCD? 

The situation regarding "exceptional configurations" has improved partly due to the finding 
[2] employing DD-HMC algorithm [3] that the numerical simulations are safe from accidental zero 
modes for large volumes. Simulations with unimproved Wilson fermions at smaller quark masses 
and lattice spacings and larger volumes have become possible with the DD-HMC algorithm. As 
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part of an on-going program [4, 5] to study the chiral properties of Wilson lattice QCD (unim- 
proved fermion and gauge actions), recently, we have demonstrated the suppression of topological 
susceptibility with decreasing quark mass in the case of unimproved Wilson fermion and gauge 
action [6, 7] where, the suppression of topological susceptibility with decreasing volume was also 
shown. In order to shed light on the mechanisms leading to these suppressions, we have further car- 
ried out a detailed study of the two-point topological charge density correlator [8]. An exploratory 
investigation of the autocorrelations of various observables with DD-HMC algorithm is presented 
in Ref. [9]. In these works we have employed ensembles of gauge configurations generated by 
means of DD-HMC [3] algorithm using unimproved Wilson fermion and Wilson gauge actions [1] 
with rif = 2 mass degenerate quark flavours. 

In this work we investigate pion mass (mj^) and decay constant, PCAC quark mass, quark 
condensate and nucleon mass in the range 290 ^ mj^ < 750 MeV. We perform qualitative chiral 
extrapolations of various observables in the range 350 ^nijc < 550 MeV. So far the simulations are 
done at two lattice spacings in the region of 0.05 - 0.07 fm where many of the modern simulations 
of LQCD with improved actions are carried out. 

2 Simulation and Observables 

Simulations have been carried out at two values of the gauge coupling correspond to j3=5.6 and 
5.8. At j3 = 5.6 the lattice volumes are 24^ x 48 and 32^ x 64 and at j8 = 5.8 the lattice volume is 
32^ X 64. The number of thermalized gauge configurations ranges from 3760 to 13646. The lattice 
parameters and simulation statistics are given in Table 1 . For all ensembles of configurations the 
average Metropolis acceptance rates range between 75 — 98%. For pion and nucleon we consider 
the following zero spatial momentum correlation functions 



where t refers to Euclidean time. For the nucleon ^i£?'2 = NN with N = {q^iCy=,qu)qu- For the 
pion &\&2 = PP\ AA\ AP^ or PA^ where P = qiYsgj (pseudoscalar density) and A corresponds to 
A4 = qiJAYsq] (fourth component of the axial vector current). Here / and j stand for flavor indices 
for the u and d quarks and for the charged pion / 7^ j. For pion we use point source and point sink 
and for nucleon we use wall source and point sink. Unless otherwise stated 20 HYP smearing steps 
with optimized smearing coefficients tti = 0.75, a2 = 0.6 and tta = 0.3 [10] are used for the gauge 
observables. 

The pion decay constant Fjt and the quark mass niq from PCAC or the axial Ward identity 
(AWI) are respectively defined, in the continuum as follows: 
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From the PP and the AP propagators 
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Table 1. Lattice parameters, simulation statistics and pion mass (nijc) in the unit of Sommer parameter ( 
ro). Here block, N2, Nc/g, T refers to DD-HMC block, step number for the force F2, number of DD-HMC 
configurations and the Molecular Dynamics trajectory length respectively. 
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Figure 1. Effective mass aiuji versus t for j3 — 5.8, k — 0.1543 and the volume 32^ x 64. 




Figure 2. Effective coefficients a^Cpp (left) and cPCap (right) versus t for j3 = 5.8, K = 0.1543 and the 
volume 32^ x 64. 

3 Computation of pion mass and decay constant, PCAC quark mass and nucleon 
mass 

The two point correlation function C{t) (with point source and point sink) may be expanded as 

C{t)=coe-'''' + cie-''''+ .... (3.1) 
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Figure 3. Effective am^^ and effective aF^^ versus t for for j3 = 5.8, k = 0.1543 and the volume 32^ x 64. 



where Mq denotes the mass of the ground state and Mi is the mass of the first excited state. The 
effective mass and the effective coefficient are given by 
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We first calculate effective pion mass {a{mji)eff) from both PP and AP correlators. Then Eq. (3.2) 
is used to determine the asymptotic value qMq of aMeff{t). Similarly we determine the asymptotic 
values of a^C^^ and a^C^^ by fitting a^C^ff and ^^C^ff respectively with the Eq. (3.3). 

The sample plots illustrating the dependence of these observables on t are given in figures 1 
and 2 for j8 = 5.8, k = 0.1543 and the volume 32^ x 64. The fitting ranges in t are 8-31 and 8-29 
for PP and AP correlators respectively. In these figures the horizontal line with error bars represent 
the asymptotic value giving the mass and coefficients. 

The effective quark mass and the effective pion decay constant are calculated by using Eq. 
(2.10) and Eq. (2.6) respectively. The sample plots illustrating the dependence of (am^^)gyy and 
{aF^^)eff on t are given in figure 3 for jS = 5.8, K = 0.1543 and the volume 32^ x 64. In these 
figures the horizontal line with error bars represent the quark mass am^^ and pion decay constant 
aF^^ calculated using the asymptotic values of the mass and coefficients from Eq. (2.10) and Eq. 
(2.6). 

We have extracted nucleon mass using wall source. Plots of the effective mass of nucleon at 



j8=5.6, lattice volume 24^ 



and K= 0.1575, 0.15775 are shown in the Fig. 4 (left) where the 



fitting ranges are t = 11-16 and 11-14 respectively. Similarly in Fig. 4 (right) we have shown the 
effective mass of nucleon at j8=5.8, lattice volume 32^ x 64 for k= 0.1543 and 0.15462 with fitting 
ranges t = 9-25 and 16-20 respectively. The effective masses are extracted from the linear fit to the 
plateau region. In table 2 we present the lattice data for aiuji, aruq, aFj^ and am^. At = 0.15466 
for j8 =5.8, the signal for effective mass of the nucleon was noisy so we do not quote a number. 
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Figure 4. Typical plots of nucleon effective mass at j3=5.6 (left) and /3=5.8 (right). 

From table 2 it is clear that finite volume effect is negligible for the pion mass and nucleon 
mass at j8=5.6 and lattice volume 24^ x 48. At /3=5.6 and lattice volume 32^ x 64 for lightest quark 
masses studied nijcL values are 4.67 and 4.23 for K = 0. 15825 and 0. 1583 respectively. Note that the 
physical volume at j8=5.6 and lattice volume 24^ x 48 is very close to the physical volume at j8=5.8 
and lattice volume 32^ x 64. Thus comparing the values of rom„ from table 1 for j8=5.8 we expect 
that finite volume effect should be negligible upto and including K = 0.15455 at lattice volume 
32^ X 64. At j3=5.8 for k = 0.15462 and 0.15466 the m^L values are 3.36 and 3.25 respectively. 
Only for the smallest quark mass (k = 0.1547) m^jL = 2.61. Hence qualitatively we expect that 
finite volume effect would be negligible for pion mass and nucleon mass for the region of parameter 
space studied except possibly for the lightest quark mass at j8=5.8. According to NLO Chiral 
Perturbation Theory (;^PT), the finite volume effect for F^^ is four times that of mj^. Thus at j3=5.8 
for K = 0.15462 and higher we can expect non negligible finite volume effect for F^^. 



4 Determination of 
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We measured the Wilson loops {W{R,T)) with temporal extents up to T = 32 (24) and spatial 
separations up to 7? = \/3 x 16 (\/3 x 12) for lattice volume 32^ x 64 (24^ x 48). 

A reasonable estimate of the static potential aV{R) is obtained by the plateau reached at large 
T of the effective potential 

{W{R,T)) 



aV^ff{R,T) = ln- 



(4.1) 



{W{R,T + l)) 

Phenomenologically, the potential V between a static quark and an antiquark at a distance r 
apart is parametrized as follows: V{r) = Vq + ar + j where a is the string tension which has the 
dimension of mass^. In lattice units, we have aV{r) = aVo + a^G ^ + af. Writing r = Ra and 
a = d/a^, we get aV{R) = aVo + a/? + |. 
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Table 2. Lattice data for am,:, arriq, aFjc and amAr. 
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Figure 5. a/ro versus aniq at j3 = 5.6 and j3 = 5.8. 



Using the expression for the perturbative lattice Coulomb potential [11, 12] 
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Figure 6. Scaling of pion mass squared (left) and nucleon mass (right). 



the parametrization of the corrected potential on the lattice reads 

^ -4l (43, 



aV{R) = aVo + a /? - ^ - 5rot 

K 



where 5rot is the coefficient of the correction term. The measured static potential is fit to the 
formula in Eq. (4.3) which corrects the lattice data for the lattice artifacts in the Coulomb potential. 
The first three terms of Eq. (4.3) now gives the continuum potential (i.e., without lattice artifacts). 
The inverse of the Sommer parameter (ro) in lattice units is calculated using 



In Fig. 4 we have plotted ^ versus lattice quark mass at j8=5.6 and 5.8. In the linear fit at 
j8=5.6, we have included k:=0.15775, 0.158, 0.158125 in the lattice volume 24^ x 48 and k:=0.158, 
0.15815, 15825, 0.1583 in the lattice volume 32^ x 64. Similarly at j8=5.8 the linear fit is done 
with k:=0.15445, 0.15455, 0.15462, 0.15466 and 0.1547 in lattice volume 32^ x 64. In the chiral 
limit, the values of ^ obtained are 0.156(4) and 0.12161(5) at j3=5.6 and 5.8 respectively. 

5 Cutoff effects 



A major source of concern with the use of unimproved Wilson fermions is the potential presence 
ff{a) cutoff effects in various observables. In this section we present the cutoff dependence of 
a variety of observables. In Fig. 6 we have plotted {r^m^Y (left figure) and r^mfj (right figure) 
versus r^Zmniq and [r^mjiY respectively using the data at j8=5.6 and 5.8. Here is the quark mass 
renormalization constant which we have calculated from the Z factors given in [13]. It is seen that 
scaling violation are negligible. 

In Fig. 7 we have plotted vqZaFji versus roZ,„mq using the data at j8=5.6 and 5.8. Data at 
j8=5.6, except at the smallest quark mass, is from lattice volume 24^ x 48 and the smallest quark 
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mass data (k"=0.1583) is from lattice volume 32^ x 64. The data at j3=5.8 is from lattice volume 
32^ X 64. Note that the physical volume at {3=5.6 and lattice volume 24^ x 48 is very close to the 
physical volume at j3=5.8 and lattice volume 32^ x 64. From the figure it appears that Fj^ data 
at j8=5.6 and 5.8 do not exhibit scaling for smaller quark mass region. Since Fj^ has larger finite 
volume effect than m^, the deviation seen at small quark mass region could be partially due to finite 
volume effects. Note that we have taken the values of non-perturbative renormalization constant 
Za at j8=5.6 and 5.8 for unimproved Wilson fermions from Ref. [13]. In this reference at j8=5.8 the 
lattice volume was 24^ x 48 and the smallest quark mass probed was greater than 75 MeV. Thus 
the systematic error in the chiral extrapolation of Za performed in the Ref. [13] at j8=5.8 may not 
be properly estimated. Hence the apparent lack of scaling exhibited by Fj^ may also be partially 
due to inaccurate determination of Za- 

In the scaling region, hadron masses are expected to be independent of the lattice scale and 
our results demonstrate that they do exhibit very small, if not negligible, cut off dependence in the 
range of lattice spacings studied. 

On the other hand the lattice renormalization constant Zy corresponding to the local vector 
current V"^'^™' which is not conserved for a / depends on the lattice spacing and is expected to 
approach unity as lattice spacing a — )• 0. In order to verify the expected lattice scale dependence of 
Zy we calculate Zy from [14] 
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Figure 8. Quark mass and scale dependence of Zy at L32T64 (left), quark mass and volume dependence of 
Zy at J3 =5.6 (right). 



Our result shown in Fig. 8 (left) exhibits the expected behaviour. On the other hand, finite volume 
effect on Zy is small though not negligible as shown in Fig. 8 (right). As already discussed, in 
order to examine the cutoff effects in pion decay constant we need accurate determination of Za in 
the chiral region which is beyond the scope of the present paper. 



6 Scale determination, chiral behaviour of nucleon mass and extraction of sigma 
term 

To extract the lattice scale, in Fig. 9 we have plotted the ratio ^ versus roniji. The data is used to 
fit the phenomenological formula 

— = airomjt - a2{rQmjtf + asirorrijt)'^ (6.1) 

motivated by baryon chiral perturbation theory. We fit the data at j3=5.8, lattice volume 32^ x 64 
and K= 0.1543, 0.15445, 0.15455 and 0.15462. At j3=5.6, we use K= 0.1575, 0.15775 in lattice 
volume 24^ x 48 and K= 0.158, 0.15825, 0.1583 at lattice volume 32^ x 64. Using the physical 
value of we get the intercept of the fitting curve and from the intercept we extract the pion 
mass in unit of ro at the physical point. Now using the value of physical pion mass and the chiral 
Umit of ^, we have computed lattice spacing a. The computed values of lattice spacings at j8 =5.6 
and 5.8 are 0.072 fm and 0.057 fm respectively. 

We have plotted dimensionfull nucleon mass versus dimensionfuU m| in Fig. 10. At j8=5.6 we 
have plotted JC=0.1575, 0.15775 at lattice volume 24^ x 48 and K= 0.158, 0.15825, 0.1583 in lattice 
volume 32^ x 64. At j8=5.8 we have plotted K= 0.1543, 0.15445, 0.15455, 0.15462 and 0.1547 in 
lattice volume 32^ x 64. To fit the data we have excluded K= 0.1575 and 0.1547 at j8= 5.6 and 5.8 
respectively. The data at K"=0.1575 is dropped because it is beyond the chiral regime and data at 
K= 0. 15475 has not negligible finite volume effect. We fit the data using baryon chiral perturbation 
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theory formula to order (m^) given in Ref. [15]: 
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For the fit we have treated Mq, ci and e\ as free parameters and set 7v= 0.086 GeV, gA= 1.256, C2= 
3.3 GeV^\ C3= -4.7 GeV^^ and A= 1 GeV. Note that the values of C2 and C3 chosen are close to their 
phenomenological values [15]. The fit along with the error is shown in Fig. 10. The values of the 
parameters Mq, ci and e\ we obtain from the fitting are Mq = 0.8 1 (4) GeV, ci = — 1 .04(5) GeV~^ 
and e[ = 1.2(1) GeV^^. Note that the ci obtained is close to the phenomenologically determined 
value [16]. 

Both from experimental and theoretical point of view, the pion-nucleon a term is of significant 
interest. A direct calculation of a involves the computation of a three point function and is beyond 
the scope of the present paper. An alternative method of calculation employs Feynman-Hellmann 
theorem utilizing the dependence of nucleon mass on the quark mass which in tern can be converted 
into the dependence on pion mass squared. With the parameters used to fit the nucleon data, we 
calculate the sigma term using the expression [15] 



2 ^8a 3 4 
-4ci m^ - ^^^^^ m^ + m^ 



^ 8C1/5 3ci 3c3 
8^1 z:^ + 



F2 8712^2 16712^2 



' 647i2MoFi ^ 4K^Fi X \ 4 4Mo) 



(6.3) 



where lj = — ( ^3 + 2 log 1 . In Fig. 1 1 , we plot a given in Eq. (6.3) together with the error 
as a function of m| and at physical point a = 0.052(4) GeV which is compatible with the currently 
available determinations [15]. 
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Figure 10. Nucleon mass versus ot| at j3=5.6 and 5.8. 
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Figure 11. Behaviour of the sigma term versus at j3=5.6 and 5.8. 



7 Chiral extrapolations of pion mass and decay constant and quark condensate 

In this section we discuss chiral extrapolation of our data obtained from the pion propagators. In 
SU(2) ;^PT at NLO [17], the quark mass (roniq) dependence of (rom^c)^ is given by 

A? 



(romjt) 



^ — Ir^rnqB 



1 



^ rln- 



167r2f2 2maB 



(7.1) 



where F is the chiral limit of the pion decay constant and B and A3 are low energy constants. 
In Fig. 12 we compare our data at j8 = 5.6 and 5.8 with a fit to the NLO ;^PT formula (solid 
Une). For clarity, the LO (dashed line) and the NLO (dot-dashed line) contributions are separately 
shown. At j8 = 5.6 and lattice volume 24^ x 48, K= 0.158 and 0.158125 and at j3 = 5.8 and lattice 
volume 32^ x 64, K= 0.1543, 0.15445, 0.15455, 0.15462 and 0.15466 are used for fitting. The 
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Figure 12. NLO chiral perturbation theory fit to {rQiriji)^ versus roZminq. 
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Figure 13. 5^ versus roZmin,, for j3 = 5.6 and i3 = 5. 



value of F = 86 MeV is taken as an input. The values of the parameters roB and roA3 obtained 
from the fit are roB = 6.20(13) and tqAs = 2.62(17). Converting these values in physical unit we 
get5 = 2651(56) MeV and A3 = 1120(73) MeV. Note that the currently available estimates of 5 
and A3 are 2112 < B < 2811 MeV and 458 < A3 < 1020 MeV [18]. Thus we find that the quark 
mass dependence of pion mass squared exhibited by our data is in accordance with NLO 
To expose the presence of the chiral logarithm, it is customary to plot 



2rnB 



1 



16712^2 2m„B 



(7.2) 



We plot 



versus r^Zmfiiq in Fig. 13. In the fit, we have used r^B = 6.2029 which is obtained 



from the fit in Fig. 12. The quark masses used in the fit are the same as in Fig. 12. 
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Figure 14. Chiral extrapolation of roF,i at j3=5.6. 



The quark mass dependence of the pion decay constant roTv in NLO ;t;PT is given by 



1 + „ o^n ln 



(7.3) 



where A4 is another low energy constant. As already mentioned, our renormalized data at 
j3 = 5 .8 has larger uncertainty coming from Za determination and also possible finite volume effects 
at the smaller quark masses. Therefore we use the data only at j8 = 5.6 for our exploratory chiral 
extrapolation. In Fig. 14 we compare our data for roFji versus r^Zmniq with NLO xPT at j3=5.6. 
Since F^ has larger finite volume effect compared to ruji, we have used data in lattice volume 
32^ X 64 in this fit. Appropriate to the chiral region, we have used only the lower quark masses 
(jc=0.158, 0.15815 and 0.1583) for the fit. Since we have fewer points to fit, we have used roA4 = 
3.1134 (A4 = 1324 MeV) (see [18]) as the input to the fit. The upper limits of the values obtained 
for F and B (see Fig. 14) are close to the currently accepted values [18]. Note however that there 
is some uncertainty arising from the value of Za- Thus in the case of it is safe to say that our 
data is not incompatible with NLO ;tPT prediction. However simulations at smaller quark masses 
with less systematic errors due to finite volume and Za determination is needed to reach a definite 
conclusion. 

Next we consider the extraction of chiral condensate from our data. In continuum limit, the 
effect of Wilson term on chiral condensate does not vanish because Wilson term is a dimension five 
operator. In other words, it renormalizes the chiral condensate additively, in addition to multiplica- 
tively. Thus in the case of Wilson fermions, a direct measurement of the quark condensate on the 
lattice and then taking the chiral continuum limit will not give the desired condensate. The well- 
known alternate procedure to calculate the chiral condensate on the lattice with Wilson fermions is 
utihzing chiral Ward-Takahashi identities [19, 20] 

(V^va) = 2mJ d^x {P"ix)P'iO)) => iWv) = 2m,^ = Cap (7.4) 
J nin 
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Figure 15. Renormalized quark condensate versus renormalized PCAC quark mass at j3 = 5.6 in unit of tq. 



using Eq. (2.10). In Fig. 15 we plot the chiral behaviour of the quark condensate (E) together with 
an NLO fit at j8 =5. 6 and lattice volume 32^ x 64. According to NLO X^T 



F^B 



1 



3Bmq 



ln(- 



, 2Bmq 



(7.5) 



The quark masses used are as in Fig. 14. In this fit we have used the chiral hmit of roF =0.20 ( 
corresponds to f = 86 MeV) as input.The value of chiral condensate obtained from the fit is Z^/^ 
= 228(8) MeV. Note that since the values of F„ have entered the determinations of condensates all 
the caveats that we have mentioned in the context of chiral extrapolation of F^^ apply in this case 
also. 



8 Summary 

We have calculated pion mass, pion decay constant, PCAC quark mass and nucleon mass in two 
flavour lattice QCD with unimproved Wilson fermion and gauge actions. Simulations are per- 
formed using DD-HMC algorithm at two lattice spacings and two volumes for several values of the 
quark mass. The cutoff effects in pion mass and nucleon mass for the explored region of parameter 
space are found to be negligible. The chiral behaviours of pion mass, pion decay constant and 
quark condensate are found to be qualitatively consistent with NLO chiral perturbation theory. 
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